Cooperative behavior of qutrits with dipole-dipole interactions 
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We have identified a class of many body problems with analytic solution beyond the mean-field 
approximation. This is the case where each body can be considered as an element of an assembly 
of interacting particles that are translationally frozen multi-level quantum systems and that do not 
change significantly their initial quantum states during the evolution. In contrast, the entangled 
collective state of the assembly experiences an appreciable change. We apply this approach to 
interacting three-level systems. 
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The description of the collective behavior of quantum 
ensembles beyond the mean-field approximation is one 
of the most challenging tasks of modern physics. The 
experiments with three-level Rydberg atoms performed 
during last decade [E H, [1] unambiguously point out the 
important role of essentially many-body phenomena in 
frozen gases with interaction among the internal degrees 
of freedom [J . Besides, the three-level systems (qutrits) 
have been considered 5] in the context of quantum infor- 
matics where, in particular, the questions of generalized 
entanglement and coherence protection Q have been 
addressed recently. Therefore, analytical results consid- 
ering dynamics of quantum assemblies 8] of three-level 
elements beyond the mean-field approach and hence al- 
lowing for the multipartite entanglement have become an 
issue of general interest. 

Here we present the exhaustive description of an im- 
portant particular class of quantum states of an assem- 
bly of N interacting identical qutrits where each single 
qutrit state remains close to the initial state with pre- 
dominantly populated middle level. Still the collective 
assembly quantum state, being essentially entangled, dif- 
fers from a product state typical of the mean-field ap- 
proximation and moreover it can considerably deviate 
from the initial state. Significant deviation of the collec- 
tive state which occurs in spite of small deviations of the 
single-particle states is known in the many-body theory 
as the orthogonality catastrophe. Our approach can be 
generalized to assemblies of arbitrary multilevel elements 
that remain close to their initial states. 

For the description we employ a technique of nilpoten- 
tials inspired by the ideas of the Glauber coherent states 
0, which has been adapted for three-level systems and 
emp loyed for the description of quantum entanglement 
[lOj . The collective state under consideration 
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can be represented in terms of the commuting operators 
tf and uf defined as the nilpotent su(3) operators 
that create the upper = jO),- and the lower | — l)j = 
uf |0)j states, respectively, by acting on the middle state 
|0)j of i-th qutrit. Here |0) = Hi denotes the initial 
state of the assembly where all qutrits are in the middle 
state fll|. The requirement (i) that each qutrit is close 



to the initial state impHes \ W\ = TtWW+ < N. The as- 
sembly state ([T]) is not normalized to unity, although the 
entanglement matrix Wij contains all information about 
the state, including normalization. It is expedient to ex- 
plicitly give the normalization factor {'^'w l^w) and the 
population ni of the upper states of qutrits 
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corresponding to the state vector eq.(IT]). It is also worth 
mentioning that the sum / — T,f^ j^^^ufWijif represents 
the tanglemeter [lOj of the state I^'vk)- Due to our ini- 
tial assumption / is not of the most generic form, i.e. it 
is lacking terms like T.f^^^^^i^uf Aijuf or Y^fj^-^tf Bijtf . 
Nevertheless / contains cross terms which ensure, ac- 
cording to the entanglement criterion [lo| , the existence 
of entanglement among the qutrits of the assembly under 
consideration. 

We now consider an assembly of qutrits with dipole- 
dipole interaction ^- Vjidjdi and find the time- 
dependent matrix W . For the dipole moment operator 
;~ -|- tf the interaction Hamiltonian reads 
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We assume that the single qutrit Hamiltonian has the 
form 



Ha = PX'^^^ + a 
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has a spectral decom- 
We also assume that 



and that the coupling matrix Vji 
position Vji = I^mQ™^™^'™- 
(3 ^ a and exclude the high frequency /S in the rotating 
wave approximation. Then the evolution operator in the 
interaction representation can be written as a functional 
integral 

exp {-iHt} = jf e-'Sl[VZ„,{t)VZ:^{t), (5) 
where the normalization constant A and the action 
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dt 



are given in terms of the components Zm{t) and Z*^(t) 
of two complex conjugated vector-function variables. 

The representation ([5]) allows one to consider dynam- 
ics of qutrits independently: each qutrit now is subject 
to an action of the single-particle time-dependent Hamil- 
tonian 



(t) = ie''^'u+ + e-'"H-)E,it) + ie-'"*u-+e'"H+)E*it) 
where 
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is an "effective electric field" depending on the functional 
variables Z„i{t). The requirement (i) justifies the em- 
ployment of the second-order time-dependent perturba- 
tion theory expression 

U,{t)^l~ij (x) dx^ J (y) J"" (x) dxdy 

(8) 

for the evolution operator of i-th qutrit initially in the 
state |0)j, which we write down in the form 



J* c'""' Ei(x)dx~itf J* c^"^ E' {x)dx 
e" r^-(y^ fc'^(y-^E:{x)dxdy ^g-j 
g- fc'^(y-^E,{x)dxdy^Qy 

It is due to this very approximation that the problem 
becomes analytically soluble. 

Substitution of ^ into (O with the allowance for the 
relation J2i V,l^^Vm^Ci„iCik SkmVm yields the action 

s = [ Zmm*^it)dt (10) 

- E /* e'"" «^™(^) + t+Z*^ix)) dx 
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which is bilinear in the field variables Z„i and Z^. This 
allows one to exactly evaluate the Gaussian functional 
integral by performing standard calculations: from the 
action S of pHl) one derives the Lagrange equations 



By substituting the solutions of these equations 

(ia sin XUm + cos XUJrn) Y.i ^m^^C™ji+ 



Zjyi{x) 



i (a + Vm) sin tujm + cos icj,, 
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(ia sin xum + ujm cos XUl^n ) Vm C, 

i (a -1- V„i) sin tuJrn + cos tuJra 



with ujm = ^Ja{a-\- 2Vm), to PH)) we find the part of 
S which depends only on the operators uf and tf but 
not on the fields Zm and Z^. The remaining part de- 
pending on the fields but not on and tf is a func- 
tional integral which gives a c-number and thus can 

be ignored. The evaluation yields exp |— i7Ji| |0) ^ 

r N ] 

cxp < J2 ^t'^ij'^t ( 1^) with the tanglemeter matrix 
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Substitution of ^ to ^ yields 
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We are now in the position to consider several particu- 
lar examples. We start with the case where all qutrits 
interact via a collective dipole moment, that is when 



V = const. In this case the matrix W has only 



one nonzero eigen value Vi 
upper state population 



NV . This results in the 
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shown in figure [T] as a function of the collective coupling 
NV and the detuning a. The maximum rate of the cre- 
ation of the qutrits in the upper states corresponds to 
a = —NV where ni — sinh^ [NVt]. This expression is 
valid as long as sinh^ [iVFt] <^ N, that is when the as- 
sumption (i) is fulfilled. One gets a deeper insight into 
the physical meaning of this result by comparing the as- 
sembly of collectively interacting qutrits with the lasing 
of an inverted two- level media [l3|. In this case the op- 
erator t+ corresponds to the photon creation operator 
af while corresponds to the excitation annihilation 
operator . As long as the total number of the emit- 
ted photons remains much smaller than the total num- 
ber of the two-level atoms these two systems are almost 
equivalent, - the main difference being that the collec- 
tive dipole-dipole interaction shifts the resonance of the 
lasing rate from a = to the point a — —NV. 

Next example is an individual dipole-dipole coupling of 
qutrits with Vji — ^^{1 — 2 cos^ 6 ji)/r^^, where r^j is the 
distance between i-th and j-th qutrits, 0ji is the angle be- 



for the extremum trajectory. tween the radius-vector and the z direction, and is 
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FIG. 1: Number ni of the qutrits in the upper state as a 
function of the coupling v = NV and the detuning a for 



the dipole moment matrix element. We consider the po- 
sition of each qutrit as an independent random variable 
and assume that the qutrits have uniform spatial den- 
sity n. The statistical properties of the random matrix 
corresponding to interaction in disordered media is 
a challenging problem which has already been addressed 
in the context of spin glasses [11 [H and cold Rydberg 
atoms [3, In particular, the distrilDution g{Vm) of the 
eigenvalues of Vij found numerically [16] has essentially 
non-analytical behavior near Vm — > 0, in contrast with 
the well-known Wigner semicircular distribution of the 
Gaussian random matrix eigenvalues. In figure [5] along 
with the results of a similar numerical work performed 
for a larger statistical ensemble and followed by a more 
accurate analytical fit to the distribution of Vm we de- 
pict the population averaged over this distribution 

g{Vm). 

Note that the numerical results with = 1 are given 
in heuristic dimensionless units a/N^/'^. Yet unknown is 
the energy parameter correctly describing the coopera- 
tive phenomena in ensembles with interaction. We 
guess that it might resemble the combination of param- 
eters jjL^ny/N where /u^n is the typical two-particle in- 
teraction while the typical parameter \/iV allows for the 
cooperative effects resembling Dicke superadience [l7| of 
two-level particles. However, unambiguous identification 
of this parameter is the subject of a more detailed future 
consideration. We also note that by assuming small de- 
viations of the qutrits from their initial state, we discard 
from the consideration the strongly interacting dimers of 
qutrits that give an important contribution to ni [2|] at 
large a. The latter effect requires a more sophisticated 
model, which would consider these dimers as the assem- 
bly elements of another type. 

The third example concerns a controlled behavior of 




FIG. 2: a) Density of the eigenvalues Vm for A'^ = 300 dipoles 
averaged over 100 random distributions in a unit cube (dots) 



and a heuristic fit g{Vm 
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green line). We set the dipole moment /i = 1. 
bution is symmetric with respect to zero, b) Number n\ of 
the qutrits in the upper state as a function of the scaled time 
tN^^'^ and the scaled detuning a/N^^'^. 



the assembly where, as earlier, the qutrits interact via 
the collective dipole, but variation of a is now possible 
during this interaction. This example shows that the be- 
havior of the three-level systems, though similar, still is 
richer than that of the two-level systems and may dis- 
play effects similar to mode beats. We consider the case 
where a, initially different from zero, can be switched 
off at t = ti, then remains zero during a time interval 
At — t2 — ti, and takes the initial value for t > t2- We 
notice the change of the symmetry associated with such a 
control: a three-level quantum system, which generically 
has the sm(3) symmetry, at a = becomes equivalent to 
a spin-1 particle possessing the su{2) symmetry. For the 
description, in the equation IjlOp for the action one has 

to replace e'"- and e'"^— ) by e'^ "(^^^^^ and e'^^ 
respectively. 

When V ^ a, the solution of the integral Langrange 
equations results in 
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sin^ [{V + a){t- At)] 
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FIG. 3: (a) Number of the particles in the upper state 
ni for time dependent a. We set a — when ti < 
t < t2 and the typical coupling V ~ y?n^/N . (b) Ratio 
ni {NV,a,t,t\,t2) /ni {NV,oi,t,ti,ti) of the numbers of the 
qutrits in the upper state as a function of the scaled time 
ti and the time difference t2 ~ t\. We set t = t2 + 0.9; 
NV = -a^ 0.3. 

for aAi ^ 1 and in 

ni ~ At^ {1 + 2 [1 - cos y (i - ta)] [1 + cos a (ti + t2)]} 

(15) 

for a At ^ 1. One can eliminate the second brackets 
in by averaging out the rapid osciUations at the 

frequency a. In figure [3^a) we depict ni (t — t2) aver- 
aged not only over these oscillations but also over the 
distribution of the collective couplings g{Vm) presented 
in figure m One can see that the population displays dy- 
ing oscillations as a function of the variable t — t2 that 
is the time elapsed after the moment t2 when the detun- 
ing a is switched back. These oscillations can be inter- 
preted as beats of the symmetric ^""'^^"'"^ and antisym- 
metric ^^"'"^^^"'"^ combinations of the upper and the lower 



states of qutrits in the interaction representation. As a 
direct consequence of su(2) symmetry, the combination 

^^^^^^"'"^ evolves while ^""'"^Z^^^ conserves when a = 0. 
The analytical solution is straightforward- though cum- 
bersome, and it does not contain any principal techni- 
cal complications. In figure EJb) we depict the ratio 
ni {NV, a, t, ti,t2) /ni {NV, a, t, ti,ti) as a function of ti 
and t2 obtained explicitly for V = const with the help of 
Mathematica package. Again one can see the population 
oscillations resulting from the quantum beats. 

We conclude by summarizing the results obtained, (i) 
The collective behavior of an assembly of three-level ele- 
ments (qutrits) admits an exhaustive analytical descrip- 
tion when each qutrit does not change significantly its 
initial quantum state. This description is based on the el- 
ements' dynamic separation with the help of a functional 
integral and invokes the time-dependent second-order 
perturbation theory followed by the exact evaluation of 
a Gaussian functional integral, (ii) In the simplest case 
where the interaction among the qutrits occurs via col- 
lective dipole-dipole interaction, the dynamics of the sys- 
tems resembles the lasing of inverted two-level media: 
for the qutrits initially in the middle state the collective 
downward transitions are accompanied by the collective 
upward transitions similar to the photon creation in a 
coherent state. The only difference is that the collec- 
tive dipole-dipole interaction shifts the resonance center 
such that the maximum transition rate occurs when the 
frequencies of the downward and upward transitions dif- 
fer by the amount of the dipole collective coupling, (iii) 
This effect persists for the regular dipole-dipole interac- 
tion and yields a two-hump frequency dependence of the 
transition rate. Positions of the maxima are found with 
the help of the collective coupling distribution, which has 
been numerically obtained for the media with 1 /r^ inter- 
action, (iv) Still three-level systems may have a more 
complex behavior than two-level ones displaying quan- 
tum beats. This effect occurs when in the course of time 
one changes the position of the middle level relative to 
the positions of the upper and the lower levels switching 
in this way between the su{2) and su(3) symmetry of 
the qutrits. 
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